Shell Stabilization of Super- and Hyperheavy Nuclei Without Magic Gaps 
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Quantum stabilization of superheavy elements is quantified in terms of the shell-correction energy. 
We compute the shell correction using self-consistent nuclear models: the non-relativistic Skyrme- 
Hartree-Fock approach and the relativistic mean-field model, for a number of parametrizations. All 
the forces applied predict a broad valley of shell stabilization around Z = 120 and N = 172-184. 
We also predict two broad regions of shell stabilization in hyperheavy elements with N ~ 258 and 
N ~ 308. Due to the large single-particle level density, shell corrections in the superheavy elements 
differ markedly from those in lighter nuclei. With increasing proton and neutron numbers, the 
regions of nuclei stabilized by shell effects become poorly localized in particle number, and the 
familiar pattern of shells separated by magic gaps is basically gone. 
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The synthesis of superheavy elements (SHE) has been 
in the focus of heavy-ion physics for more than three 
decades. The last few years have seen significant progress 
in our quest for reaching the region of long-lived SHE. 
Light isotopes of the elements Z = 110-112 have been 
safely established at GSI Darmstadt and JINR Dubna 
These isotopes are expected to be strongly de- 
formed thanks to the (predicted) deformed shells Z = 108 
and N — 162 (see Refs. and references quoted 

therein). These new nuclides could be unambiguously 
identified by their characteristic cv-decay chains lead- 
ing to already known isotopes. Even heavier and more 
neutron-rich nuclides have been announced just recently 
by GSI (5|gll0) §; Dubna (?f?112, ^1^114, and 
??§116) §; and Berkeley (f% 118) [|. The a-decay chains 
of those nuclei cannot be linked to any known nuclides 
as they end with fissioningnuclci. However, these results 
yet need to be confirmed ||Jl^,|ll| . 

The mere existence of SHE relies on quantum mechan- 
ics. According to the classical liquid-drop picture, all su- 
perheavy nuclei should be unstable against spontaneous 
fission — due to the huge Coulomb repulsion. How- 
ever, additional stabilization of binding energy is possi- 
ble thanks to shell effects which generate local minima in 
the nuclear potential energy surface in the regions where 
the level density around the Fermi level is lowered. The 
detailed energy balance between the local minima is dic- 
tated by the distribution of spherical single-particle or- 
bitals. In some cases the minima are sufficiently deep to 
stabilize the nucleus against spontaneous fission; the de- 
lay in the spontaneous fission half-lives due to the shell 
effects can be as much as 15 orders of magnitude for 
Z > 106 @. 

The half-lives of the known isotopes of elements with 
Z > 105 are predominantly limited by a decay and de- 
crease from 0.9 s for ff|l06 to 0.2 ms for ^112. These 
isotopes decay mostly by groups of successive a particles. 



Although shell corrections strongly influence Q a values, 
there is no simple correlation between the magnitude of 
shell effects and a-decay half-lives For instance, if 

the shell corrections are nearly constant in a broad region 
of particle numbers, this will have very little influence on 
Q a values; hence on T a . 

The subject of the present paper is the shell stabiliza- 
tion quantified in terms of the shell correction -E s hcii- It 
is obtained from decomposing the self-consistent binding 
energy E tot as 



E t , 



E + E, 



shell i 
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where E is the average energy that changes smoothly 
with particle number. (In microscopic-macroscopic ap- 
proaches, E is approximated by the liquid drop or droplet 
model energy.) The shell stabilization of SHE as pre- 
dicted by macroscopic-microscopic models has already 
been extensively discussed in the literature ^,Q^5). It 
is the aim of the present study to analyze shell correc- 
tions in the superheavy region in the framework of self- 
consistent mean-field models. To this end, we apply the 
two most widely used approaches, namely the Skyrme- 
Hartree-Fock (SHF) theory and the relativistic mean- 
field (RMF) theory. (For a brief overview, see Ref. Jl6||.) 
In the previous paper by Kruppa et al. |L7| , the shell en- 
ergy extracted from self-consistent single-particle spectra 
along a few selected isotopic and isotonic chains was dis- 
cussed. Here we present a large-scale survey of spherical 
shell energies throughout the whole landscape of conceiv- 
able SHE. 

There is a world of different parametrizations for SHF 
as well as RMF. They agree more or less in their perfor- 
mance for stable nuclei but can yield differing predictions 
in extrapolations. For example, magic shell closures in 
SHE are at variance 1 13 17|-[l9|]. Interestingly, it has been 
concluded in Ref. that both the SHF and RMF calcu- 
lations were internally consistent. That is, all the Skyrme 
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FIG. 1. Total shell correction (sum of proton and neutron shell corrections; in MeV) calculated for spherical even-even 
nuclei. The thick solid lines denote two-particle drip lines. Black squares mark nuclei calculated to be stable with respect to /3 
decay. White color indicates nuclei with positive shell corrections, black color denotes nuclei with E a he\\ beyond —12 MeV. 



models with conventional spin-orbit force predicted the 
strongest spherical shell effect at N — 184 and Z = 124, 
126, while all the RMF forces clearly preferred N = 172 
and Z — 120, and SHF parameterizations with relativis- 
tically extended spin-orbit interaction were in between. 

In view of these differences, we consider several 
parametrizations per model. We use a selection of forces 
which have been found earlier to represent the whole 
range of possible predictions for spherical shell closures. 
For the RMF, we consider the parametrizations NL3 pctl 
and NL-Z2 ]l9j . Both represent recent fits which perform 
very well with respect to global ground-state properties 
but differ in detail. NL3 reproduces well isotopic trends, 
while NL-Z2 also fits the electromagnetic nuclear form 
factor. For the SHF, we consider SkP |21| as a represen- 
tative for a force with effective nucleon mass m* /m — 1, 
leading to a comparatively large density of single-particle 
levels. All other SHF forces employed here have smaller 
effective masses around m* /m « 0.7. SLy6 j22| was ad- 
justed with particular emphasis on isotopic trends and 
neutron matter. SkI3 and SkI4 [||| employ an extended 
form of the spin-orbit force which was found to be nec- 
essary for a description of isotope shifts in heavy Pb iso- 
topes. For SkI3, the extension was restricted to map the 
spin-orbit structure of the RMF as closely as possible. 

The shell energies are computed using the same pre- 
scription as outlined in Ref. |17|| . This procedure, based 
on the Green's function approach to the level density, 
is better suited for the calculation of weakly bound sys- 
tems than the traditional approach. This is important in 
the context of nuclei considered here since some of the 
predicted regions of shell stability lie very close to the 
proton drip line. In our calculations, we include a large 
space of single-particle states up to 50 MeV above the 
Fermi energy. Since most of these states are continuum 
(positive-energy) states, the contribution from a parti- 



cle gas (treated in the same numerical box) has to be 
removed |lj] . Pairing correlations are ignored. 

The calculations are restricted to spherical symmetry. 
Consequently, the calculated shell corrections represents 
in most cases an upper bound. In many cases, deforma- 
tion does provide an additional binding, i.e., pulls to even 
stronger shell stabilization. Unfortunately, large-scale 
symmetry-unconstrained calculations of shell effects in 
the SHE are currently beyond our reach. This is because 
triaxial Jl3| and reflection-asymmetric j24| shapes must 
be considered together with more exotic topologies (e.g., 
bubble (25) , toroidal, and rod structures) which might be- 
come favored for the heaviest systems investigated here. 

For nuclei up to Z = 82, -E s heii is always sharply peaked 
at shell closures 26 2T|] . Figure |l| shows an example from 
the Sn region. The magic numbers Z = 50, N — 50, and 
N = 82 clearly stick out for all forces. One basically sees 
sharp stripes along magic proton and neutron numbers. 
All forces show the same magic numbers, but there are 
some differences in detail. The overall strength of the 
shell effect seems to scale with effective mass, with SkP 
giving the smallest shell energies, while they are most 
pronounced for NL3. 

The systematics of shell corrections change dramati- 
cally when going to SHE; see Fig. ||. Instead of narrow 
stripes of large -E s h e ii localized around magic numbers, 
all forces employed predict a wide area of shell stabiliza- 
tion which spreads over all shell closures predicted by the 
various forces. Unlike in normal nuclei, large shell correc- 
tions in SHE can appear slightly away from shell closures. 
This was, in fact, already recognized in macroscopic- 
microscopic models p jjljjisf . As a consequence, the sig- 
nificant differences seen in the prediction of magic shells 
through various binding-energy indicators (such as 82 q 
Jl8]|) are much mellowed by the generally softer pattern 
of the shell energy. We will discuss that aspect in more 



2 




160 180 200 

n Neutron Number N 

FIG. 2. Same as Fig. hi except for superheavy nuclei around the expected island of stability around Z — 120 and N = 180. 
The scale for E s h E u is the same as in Fig. 0. See Ref. @ for the individual shell correction of protons and neutrons along cuts 
through these maps. 



detail below. One of the common features seen in Fig. g 
is that the region of nuclei with largest shell corrections 
forms a triangle with the base at N — 184 and the outer 
corner at f y| 120. This happens because the existence of 
the N = 172 neutron subshell is strongly coupled to the 
proton subshell closure at Z — 120 JTgj ] . 

At second glance, however, one also sees differences 
among various parameterizations. More stable nuclei 
are found above Z — 120 for SkP which predicts the 
proton gap at Z = 126, while the center of gravity is 
clearly shifted below Z = 120 for SkI4 with its strong 
Z = 114 shell effect. The other forces reside in between. 
A somewhat different bias is also seen for the extension 
in neutron direction. SkP predict strong shell effects 
for a number of nuclei with TV > 184, while other forces 
fill basically the landscape between the two magic num- 
bers N — 172 and N = 184. There are also some differ- 
ences concerning the overall area of the stabilized region. 
For instance, NL3 makes it much smaller than all other 
forces. Of course, for a quantitative discussion, one needs 
to account for deformation effects which will serve to ex- 
tend the island of shell stabilization. For example, the 
well-known region of deformed shell-stabilized SHE lo- 
cated around fg§Hsios is missing in Fig.0, as well 
as the deformed shell closure at N — 174 |28 29]. 

Figure ^ takes a daring glance to even heavier nuclei 
with Z > 126 (hyper heavy elements). There are several 
broad valleys of spherical shell stability showing up that 
extend around the actual shell closures. This is a com- 
mon feature of very heavy nuclear systems. Differences 
between the forces, however, grow dramatically in the hy- 
perheavy region. The upper limit for the plots is chosen 
to stay below the region where the novel topologies, i.e., 
semi- bubbles and bubbles, are believed to coexist [ B5[ . 
As the RMF predicts the lower border of this transitional 
region at smaller values of Z than SHF, we set different 
upper limits of the displayed area in both approaches. By 



inspecting Fig. |^ carefully, one can see rather large differ- 
ences in -E s heii predicted in different calculations. That is 
no surprise because fine details of shell structure play an 
increasing role with increasing nuclear size. Nevertheless, 
there still remains an overall agreement concerning the 
position of the regions of stability, around N = 258 and 
around TV = 308. The RMF parametrizations are more 
pessimistically predicting only faint effects, while SHF 
produces a strong stabilization at N = 308. 

It is an open question whether for these hyperheavy 
elements the actual shell corrections are sufficient to pre- 
vent (or significantly slow down) spontaneous fission. 
One would expect that, with increasing Z 7 the Coulomb 
force would act increasingly against stability. In particu- 
lar, the role of triaxial or reflection asymmetric degrees of 
freedom must be considered when assessing the stability 
of hyperheavy nuclei to fission. 

Figure |I] shows the single-particle spectra for three typ- 
ical nuclei from the three regions of large shell correction 
discussed in this paper. For 132 Sn, proton and neutron 
magic gaps appear in all models. The patterns of single- 
particle levels are significantly different for the two re- 
gions of SHE. Firstly, with increased mass, the overall 
level density grows as oc A 1 / 3 . Secondly, no pronounced 
and uniquely preferred energy gaps appear in the spec- 
trum. This shows that shell closures which are to be 
associated with the gaps in the spectrum are not robust 
in that region. Tiny changes in, e.g., spin-orbit prop- 
erties can shift the gaps substantially; see (l^Jl^l for a 
thorough discussion. Protons are more sensitive in that 
respect than neutrons. 

Interestingly, similar problems are encountered in 
atomic calculations of the electron shell structure of SHE 
pO| . Due to the large density of valence electron shells, it 
is extremely difficult to make robust predictions of chem- 
ical properties of SHE. 

A close inspection of Fig. allows for a rather good 
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and H but for larger neutron and proton numbers (up to Z — 200 for SHF and Z — 170 for 
RMF) calculated with the subset of forces as indicated. The energy scale is kept the same as in the other figures. 



understanding of the shell-correction pattern discussed 
above. For instance, in ^126 there appear low-j single- 
particle orbitals at the Fermi surface (3pi/ 2 and 3p 3 / 2 in 
the protons between Z = 120 and 126 and isi/2, 3d 3 / 2 , 
and 3d 5 /2 in the neutrons between N = 172 and 184). 
The low (2j+l) degeneracy of these shells gives rise to re- 
duced single-particle level densities; hence to a large neg- 
ative shell-correction energy for a whole range of neigh- 
boring nuclei. 

For hyperheavy nuclei, the level density is even higher, 
and no strong shell closures are predicted in most models. 
It is only in SLy6 and SkI3 that a relative large N = 308 
gap is predicted, bounded by high-j shells. This is con- 
sistent with Fig. |jj which shows a rather strong neutron 
shell effect at N = 308 for these two forces. 



In summary, we have investigated the spherical shell- 
stability of superheavy elements using state-of-the-art 
self-consistent models. The investigation of the system- 
atics of shell energy reveals a new feature when going to 
very heavy systems. Beyond Z = 82 and N — 126, the 
familiar localization of the shell effect at magic numbers 
is basically gone. Instead, the theory predicts fairly wide 
areas of large shell stabilization. Consequently, there is a 
good chance to reach shell-stabilized SHE experimentally 
using a range of beam-target combinations. 

The disappearance of a familiar pattern of magic num- 
bers and the appearance of broad valleys of shell stability 
is due to (i) the rather large single-particle level density, 
and (ii) the appearance of many low-j shells around the 
Fermi level. This explains the robustness of the shell cor- 
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FIG. 4. Single-particle spectra of protons (top) and neutrons (bottom) for 132 Sn, ^§4 126 and 30! 164. The single-particle 
energies are taken relative to the Fermi energy predicted by SLy6. The line thickness of each level is proportional to the 2j + 1 
degeneracy of the state. The inset in the rightmost panel shows the predicted radial neutron distributions. No bubble structure 
is predicted for 30! 164. 



rcction in a rather large range of SHE and at the same 
time the volatility of magic shell closures. 

The results presented here have to be taken with a 
grain of salt. Large shell correction is a neccesary, but 
not sufficient, condition for the appearance of long-lived 
SHE. The calculations presented in this paper serve as 
a starting point for subsequent studies of deformed shell 
effects and fission barriers in superheavy and hyperheavy 
nuclei. 
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